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I. INTRODUCTION 



An issue of the tunneling time, i.e., the time associated with the passage of a particle 
through a tunneling barrier, has been discussed in many theoretical studies HH-ITf and is 
not settled yet. This difficulty arises mainly from the fact that time is not an observable 
represented by a self-adjoint operator but just a parameter in quantum mechanics. 

In our previous paper [OT, we proposed a new method to evaluate the tunneling time, 
using Nelson's approach of quantum mechanics [0. Our aim then was to treat tunneling 
effects in a detailed time-dependent and fully quantum mechanical way, as any theoretical 
expression of the tunneling time must be tested by experiments which are feasible at present 
and in near future. 



As discussed in the reference p0| , the Nelson's approach of quantum mechanics has 
several advantages to study the tunneling time, a few of which are listed below. 

First of all, this approach using the real time stochastic process enables us to describe 
individual experimental runs of a quantum system in terminology of "analogue" of clas- 
sical mechanics. This is true even in the tunnel region where classical path is forbidden. 
From sample paths generated by the stochastic process we obtain information on the time 
parameter, in particular, the tunneling time. 

As a matter of course, the whole ensemble of sample paths gives us all the same results 
as quantum mechanics in ordinary approach does, e.g., expectation values of observable, 
transmission and reflection probabilities in scattering problem and so on. It is important 
for us to note that in scattering phenomena (those without bound states) the transmission 
and reflection ensembles are defined unambiguously, namely each sample path is classified 
distinctively into either transmission ensemble or reflection one. 

We need to accumulate a sufficient number of sample paths in numerical simulations. In 
thick or/ and high potential case the transmission probability is low and consequently we have 
a difficulty that a number of sample paths belonging to the transmission ensemble is also low 
when each sample path is followed in forward time direction. However, in Nelson's approach 
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there is not only the forward Langevin equation but also the backward Langevin equation 
(see (Q) below), both being equivalent to each other in physical results. The difficulty above 
is avoided when the backward Langevin equation is employed. 

Taking account of these advantages, we have developed a theoretical formulation of time- 
dependent description of tunneling phenomena based on the Nelson's stochastic approach in 
j20fl . Numerical simulations for a one- dimensional square well potential barrier model were 



demonstrated. An important result about the tunneling time then is that there are the three 
characteristic times, i.e., the passing time and the hesitating time, and their sum, the inter- 
acting time. The probability distribution of these three times were calculated numerically. 

Our previous study treated only a quantum system of a single particle under a simple 
potential. But realistic experimental situations are much complicated. Naturally we are 
tempted to extend our previous formulation to more general scattering phenomena. In 
this paper we consider cases in which transition processes into other channels or absorptive 
processes takes place during scattering processes, and look into these effects on the tunneling 
time. 

Processes of transition into other channels and absorption are described by channel 
coupling and optical potential (complex potential), respectively, in the ordinary quantum 
mechanics using the Schrodinger equation. So far it is known that the Nelson's formulation 
is equivalent to the Schrodinger equation only for a one-body problem with a single channel 
and a real potential. The purpose of this paper is to generalize the Nelson's stochastic quan- 
tization so that it can deal with multi-channel coupling and/or optical potential problems. 
As will be shown below, one can construct such generalized formulations of the Nelson's ap- 
proach with additional stochastic jumping processes. These theoretical formulations allow 



us to perform numerical simulations of stochastic processes as before ||20|| . This way we can 
investigate the effects of transition into other channels or absorption on the tunneling time. 

The paper is organized as follows: In the next section the original Nelson's quantum 
mechanics is reviewed briefly for later relevance. We propose a formulation of the Nelson's 
approach, generalized to a quantum system with channeling-coupling in Sec. 3. The formu- 



lation of Sec. 3 hints us how to develop a formulation for optical potential, which is shown 
in Sec. 4. In Sec. 5 numerical simulation for square well potential model, using the formu- 
lations in Sections 3 and 4, are demonstrated, and physical implications of these results are 
analyzed. The final section is devoted to summary and some comments. 



II. BRIEF REVIEW OF NELSON'S QUANTUM MECHANICS 

We start with a brief review of the original Nelson's quantum mechanics which consists 
of two basic conditions, i.e., the kinematical condition and the dynamical one. 

The kinematical condition is given by the Ito-type stochastic differential equation: There 
are two ways to express it, depending on forward or backward time direction. Explicitly we 
have for forward time evolution, 

dx(t) = b(x(t),t)dt + dw{t), (1) 

and for backward time evolution, 

dx(t) = b*(x(t),t)dt + dw*(t). (2) 

The dw(t) is the Gaussian white noise (representing the quantum fluctuation) with the 
statistical properties of 

h 

< dw(t) >= and < dw(t)dw(t) >= —dt, (3) 

m 

and the same properties for dw*(t) as in ([3|). Here < • • • > means a sample average. It is easy 
to show that for these two Langevin equations hold the following Fokker-Planck equations 
for the distribution function P(x,t) of the random variables x(t), 



dP{x,t) 

dt 
dP(x,t) 

dt 



d , . , h d 

-—b(x,t) + 
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dx ' 2m dx 

d , , . h d 2 
—b*{x,t) + 



P(x,t) (forward in t), (4) 
P(x,t) (backward in t). (5) 



dx 2m dx 2 

Thus a pair of equations (JTJ) and (Q) is mathematically equivalent to a pair of equations 
and (13). We get an osmotic velocity, u from the sum of Eqs. (Eh and (M) as 



_ b-K _ 

U ~ 2 ~ 2mPdx ^ 
under the boundary condition of 

P(x^oo,t)^0. (7) 



The subtraction (H) from (Eh gives 



<9P d 



where v is a current velocity 



v 



0) 



The elimination of P(x,t) from (H) and (§) leads to an equation called the kinematical 
equation, 

du h d 2 v d 

m=-2^d^-di- {uv) - (10) 

The dynamical condition is expressed through the "mean time derivatives" introduced 
as follows: The "mean forward time derivative" Df(t) is defined as 

Df(t) ee A hm o < f(t + A £~ m \m (s < t) fixed >, (11) 

and the "mean backward time derivative" _D*/(t) is defined as 

DJ(t) = lim < ~{'*~ At \ f(s) (s > t) fixed > . (12) 

The "mean balanced acceleration" is introduced through the definitions of (|TT|) and (|l|) as 

a(x(t),t) = — — —x{t). (13) 

Note that this definition can be rewritten as 

h d 2 u 1 d , 9 9n dv 
afet) = -^ + 2^ ( " (14) 
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from (|l|) and (0) with (]6|) and (|9]). The dynamical condition is nothing but the classical 
Newton equation to this "mean balanced acceleration" a(x(t),t), that is, 

dV 

ma(x,t) = -—, (15) 

from which we derive the "Newton-Nelson equation" , 

dv h d 2 u dv du 1 dV 
dt 2m dx 2 dx dx m dx 

because of (|l4l) . 

Summarize the mathematical structure of Nelson's quantum mechanics. The two basic 
equations, (|II]) from the kinematical condition, and ([il|) from the dynamical condition, form 
a set of simultaneous equations for two unknown functions u(x, t) and v (x, t), or equivalently 
b(x, t) and t). Then we can determined the ensemble of sample paths or the distribution 
function P(x,t). Although it is practically very difficult to solve these equations directly 
due to their nonlinearity, one can easily show the equivalence between this approach and 
the ordinary approach of the Schrodinger equation. This fact tells us that one can solve the 
problem by means of the wave function much more easily. The equation 



d_ 

dx 



.h 1 dip' 1 fh\ 2 1 <9V _ 1 v 
mip' dt 2 \ m J ip' dx 2 m 



(17) 



follows from the combination of (|i0|) +i (jig) , where 



h 1 dip' 

u + lv= r (18) 

m ip dx 

Equation (|P7|) clearly shows the relationship between ip' and the wave function ip as the 
solution of Schrodinger equation 

.dip ( h 2 d 2 \ , 

namely 

ip{x,t) =ip'(x,t)exp{ — — / r}(s)ds) (20) 



with an arbitrary function of t, i](t), which has no physical relevance. It is easily seen from 
this proof of the equivalence that one has the expressions for b(x,t), b*(x,t) and P(x,t) in 
terms of ip(x, t), 

h d 

b(x,t) = — (Im + Re)— lnip(x,t), (21) 
m ox 

% d 

bJx,t) = — (Im-Re) — \nip(x,t), (22) 
m ox 

P(x,t) = \^(x,t)\ 2 . (23) 



III. STOCHASTIC FORMULATION FOR QUANTUM SYSTEM WITH 

CHANNELING COUPLING 



We now turn to generalize the above Nelson's approach to a system with a channel 
coupling. For simplicity, consider the 2-channel Schrodinger equations = {1,2}) 



. d . . . 

in—ipi{x,t) 

ot 



K_cP_ 

2mj dx 2 



+ Vu(x,t) ipi(x,t) + V ij (x,t)ij; j (x,t), 



(24) 



with 



V ■ = V* 



(25) 



Here and below the dummy index does not imply taking a sum. As will be seen, the 
generalization of the formulation in this section to N-channel case (N> 2) is straightforward. 

Consider the Fokker-Planck equations in the stochastic formulation, corresponding to 
(E3). First we require a natural extension of (El) to the present case, 



Pi(x,t) = Mx,t)\ 2 . 

The diagonal parts (the kinetic energy and Vu terms) in 
as in the previous section. The Schrodinger equations 
suggest the following equations for Pi(x, t), 



(26) 

are expected to be dealt with 
and their complex conjugates 



Ot 



d_ b + h d 2 
dx ' 2m,- dx 2 



>,t) 



Pi(x,t) (forward in time), (27) 



dPj(x,t) 
dt 



^ (x ' t) + 2^^ + ^ (x ' t) 



Pi(x,t) (backward in time), (28) 



as Pi(x, t) increases or decreases, due to the potential Vij causing transitions between i and 
j, at the rate of the absolute value of 

2, 



(29) 



Although the sum of ( 27|) and (28) leads to @ with the index i, 

bi — b*i hi dPi 



Ivdi Pi dx 



(30) 



their difference provides us with 



instead of (l8f), where 



bi + 



(31) 



(32) 



As a result, eliminating Pi(x,t) from ( |3~0D and (plj), one derives the following kinematical 
equation 

h d 2 Vj d , 



duj 
dt 



2m,- ox 



(33) 



2rrii dx 2 dx 
instead of (|T0|) . 

Here arises a natural question what are stochastic differential equations corresponding 
to the Fokker-Planck equations in (p7|) and (p8|) , just as (JT|) and @ correspond to (|J) and 
(§). Apparently we need two random variables Xi(t) (i = 1,2), which are assumed to be 
subject to the stochastic differential equations, similar to ([j]) and (fj), 



d,Xi(t) = bi(xi(t),t)dt + dwi(t) (forward in time), 
dxi(t) = b*i(xi(t), t)dt + dw*i(t) (backward in time), 



(34) 
(35) 



with the properties for dwi(t) and dw^(t), 



h 

< dwi(t) >= 0, < dwi(t)dwj(t) >= — Sijdt 

< dw*i(t) >= 0, < dw*i(t)dw*j(t) >= — 5ijdt. (36) 

TYli 

As is easily seen, a naive interpretation of these independent stochastic differential equations 
leads only to the Fokker-Planck equations in ( |27| ) and fl28|) without the terms proportional 
to Wu-*j). An additional mechanism to take account of the quantum jump between i and 
j represented by the terms involving Wu^j\ is necessary. For this purpose we supplement 
([34]) and ([35]) with a stochastic jumping process between i and j. Thus we attempt below 
the formulation for the two random variables Xi(t), subject to the stochastic differential 
equations (|3~4"D and (|35| ) combined with a stochastic jumping process in the following way. 

The "dynamical" rule to determine how each sample path Xi(t) changes its index (i = 
1 —>■ 2 or vice versa) during passage of time is described by the following random jumping 
process (Fig. 1): At each time a dice is cast, independently of the stochastic equation ( |34"D 
and (|35|), and each sample path either keeps or changes its index at a certain rate. For 
forward time direction, we have the rule in case of W^j) > (i ^ j), 



Xi(t) 



Xj(t + dt) with the probability of W(i-*j)(xi(t), t)dt, 
Xiit + dt) with the probability of 1 — W^j){xi{t) , t)dt 
x j{t) — > x j(t + dt) with the probability of 1, (37) 



and the rule in case of W^j\ < 



Xj(t) 



Xi(t + dt) with the probability of — W^j)(xj(t), t)dt, 

Xj(t + dt) with the probability of 1 + W^j)(xj(t), t)dt 
Xi(t) — > Xi{t + dt) with the probability of f . (38) 



Likewise, the rules for backward time direction state that in case of W^j) > 



Xj(t) 



Xi{t — dt) with the probability of W^j)(xj(t), t)dt, 

Xjit — dt) with the probability of 1 — Wa^j\(xj(t),t)dt 
Xi(t) — > Xiit — dt) with the probability of 1, (39) 



and in case of W^j) < 



Xi{t) 




Xj(t 

Xi(t 
Xj(t 



dt) with the probability of — W^j)(xi(t), t)dt, 

dt) with the probability of 1 + W^j){xi{t), t)dt 
dt) with the probability of 1. 



Xj(t) 



(40) 



According to the rules of the random jumping process above, the behavior of each sample 
path is illustrated as follows: For forward time direction, a sample path starts from Xi(ti), 
develops according to fl34[) with i for a while, and when a chance comes, it changes its index 



process from X{ to Xj is allowed and the reverse process is forbidden where Wu^j\ > 0, and 
vice versa where W^j) < 0. The jumping processes may be repeated or may not occur, 
depending on the sign and magnitude of W(i->j). Sample paths show similar behavior for 
backward time direction. 

It is remarked that Xi(t) is generally a functional of both of dwi(s) and dw2(s) (s < t) 
(or dw*i(s) or dw^(s) (s > t)) as it may repeat jumps between i = 1 and % = 2 in the past 
(in the future). Due to changes in the index for each sample path, there are several types 
of averages which are distinguished from each other carefully. It is convenient to introduce 
notations for conditional averages. The simple average < • • • > should be taken over both of 
dwi(s) and dw 2 (s) (s < t). To represent a physical average of the i-state at t, we introduce 
a notation of 



where the average on the left-handed side implies a conditional average only over sample 
paths, labeled by i at t. This average should be expressed in terms of the probability 
distribution Pi(x,t) as 



from i to j and follows ( |34"D with j until the next jumping process takes place. The jumping 



< f(x(t)) > {xi{t)} =< f(Xi(t)) > 



(41) 




(42) 



The notation f(x(t)) ^ > { Xl (t)}u{x 2 (t)} owns its trivial interpretations, 
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< f(x(t)) ^>{ Xl (t)}U{x2(t)} = < f(x(t)) > . 



(43) 



Furthermore, conditional averages with different times such as <C f(x(t)) ^>{ Xi (t+dt)}n{xj(t)} 
can be introduced: This example represents the average only over the sample paths which 
have the index j at t and i at t + dt. 

Let us now evaluate the time derivative of the physical average <C f(x(t)) ^>{ Xi (t)}- For 
forward time direction, using appropriate conditional averages, we write down as 



d 
dt 



< f(x(t)) > {XiW} 



1 

~dt 
1 

~dt 



< f(x(t + dt)) ^> {Xl (t+dt)} - < f(x(t)) 

< + dt)) - f(x(t)) »{ Xi (t+dt)}n{x i (t)} 
+ < + dt)) »{x i (t+dt)}n{* i (t)} 

- < /(a:(t)) »{ a!j (t+«ft)}n{a! i (t)} 



(44) 



The three terms here are manipulated as 



< /(x(t + dt)) - >{x i (t+dt)}n{x i (t)} 



df(x) 



dx 

dfjx] 
dx 



dx(t) + 

a;=a;(£) 

bi(x(t),t)dt + 



Id 2 f{x) 



2 drr 2 
d 2 /(x) & 



{dx{t)) 2 + o(dt 3/2 ) > {a;i (t+^)}n{x l (t)} 



<ir 2 2m 



x=x(t) 

dt ^>{ Xi {t+dt)}n{ Xi (t)} 



df(x) , , , N . , d 2 f(x) k , . , 9 . 

=<< dx^ bMt) ' t)dt+ ^x^2m~ dt >>{ ^ {t)} + ° (rft } 



= cit / <ir 



9 . , , ft <9 2 x 



y dxf(x) l-—bi(x,t) 



+ 



2m,- <9x 2 



Pi(x,t) + o(dt 2 ), 



(45) 



< + dt)) ^> {Xl (t + dt)}n{ Xj (t)}= -dt J dx f(x)W {i ^ j) (x,t)P i (x,t)9(-W {i ^ j) (x,t)) + o(dt 2 ), 



(46) 



and 



< f(x(t)) > { ^( t+dt ) } n{x i (t)}= * y da;/(a;)^ (l ^ ) (x,t)P l (a;,t)^(^ j) (x,t)) + o(dt 2 ), (47) 
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respectively, from (@), ([36]), (|37D, (|38|) and (f4|). Collecting fl4|)~(||7|), we obtain the correct 
time evolution of (0). This shows the equivalence between ( |27| ) and the stochastic equation 
(0) supplemented with the stochastic jumping process (|37| ) and (|38|). Likewise one can 
show the equivalence between (|28| ) and the stochastic equation ([35]) supplemented with the 
stochastic jumping process ( PP| ) and (|I]). 

We need some careful treatment on the dynamical condition in the present case. For the 
equivalence between the Nelson's approach and the Schrodinger approach, the dynamical 
condition is desired to have the form of 

1 dVu 



dvi h d 2 Ui dvi dui 
dt Ira,-, dx 2 1 dx 1 dx 



(48) 



rrii dx 

Here we introduce a "quantum potential" Vu which is to include the effect of channel coupling 
as well as the usual potential Va. The simplest way to achieve this equation is to define the 
"mean balanced acceleration" a« through the "mean (forward and backward) time derivative" 
as usual but for the stochastic process without any jumping process. We simply consider a 
stochastic process governed by (p4|) all the time, and denote Xi instead of Xi to distinguish 
them from each other. There is no mixing of dwi and dwj in Xi, contrary to Xi. For each 
Xi(t) we define the "mean balanced acceleration" aj(Xj(t),t), and the "Newton" equations, 

8% 



miai{Xi(t),t) 



becomes 

The combination of the equations 



)+*( 



dXi 



derives 



(49) 



d_ 

dx 



m ipl dt 2 I mi 



1 d 2 ^j 
ifj'i dx 2 



-1 {% - %w t 

rrii 



(50) 



where the relation, 



Ui + ivj 



h l di/j'j 

rrii ipl dx ' 



(51) 



is used. If we shift the function ■?/>• to 



ipi{x,t) = ip'i{x,t) exp(- 



zrrii 



r](s)ds), 



(52) 
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choose the "quantum potential" as 

V u = V u + Re 

and use the relation 



m 2 



(53) 



(54) 



we can reproduce the Schrodinger equations fl24]) . By the use of the flSll) and (|52|), the 



relations, 



h d 
bAx,t) = — (Im + Re)— — InibAx, t) 
rrii ox 

h d 
Ki(x,t) = — (Im-Re) — hi ipi(x,t) 
rrii ox 



(55) 
(56) 



and (|26|), are established again. 



IV. STOCHASTIC FORMULATION FOR QUANTUM SYSTEM OF OPTICAL 

POTENTIAL 



In this section, let us formulate the Nelson's stochastic approach to a system of a single 
degree of freedom described by an optical potential. Then the Schrodinger equation with an 
imaginary part of potential, denoted by ill (a physically relevant situation, i.e., an absorptive 
process corresponds to U < 0), is written down as 



ih 



dip(x, t) 



h 2 d 2 



+ V(x,t)+iU(x,t) ip(x,t). 



(57) 



dt \ 2m dx 2 

The formulation in the previous section suggests a method to establish a stochastic 
formulation for this Schrodinger equation. The analogy between the channel coupling model 
and the present model becomes apparent when we attempt the Fokker-Planck equation 
corresponding to fl57D in the form of 

dP(x,t) 



dt 
dP(x, t) 
dt 



d h d 2 2U 
b H 1 

dx 2m dx 2 H 

d h d 2 2£T 
dx * 2m dx 2 Ti 



P(x,t) (forward in i), 
P(x,t) (backward in t). 



(58) 
(59) 
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Equations (|58| ) and (^) are compared with ( p7| ) and (p8]), then both are quite similar to 
each other with the correspondence between 2U /Ti and —Wu^jy 

While the sum of (|58|) and (|59|) is given by (0), their difference leads to 



instead of (|8|). From (|5|) and (|5"UD, follows the kinematical equation 

du h d 2 v d . , ^ 1 9 ^ 
<9t 2m <9x 2 dx mdx , 

instead of fllOl). 

The additional term in (|58| ) simply describes production (absorption) effects for U > 
([/ < 0), which one may put in such a way that the production (absorption) process is a 
transition from a "unphysical" sector to a "physical" one (from a "physical" sector to a 
"unphysical" one). At this point the analogy between the previous section and this section 
is helpful to find stochastic processes equivalent to the Fokker-Planck equations in (|S8|) and 
(p)5|) : We consider the two random variables x p (t) and x u (t) for "physical" and "unphysical" 
sectors, respectively, and stochastic jumping between them occurs according to certain rules, 
which will be specified below. In contrast with the channel coupling case with the index i, 
the stochastic differential equations for both of x p (t) and x u (t) can be common. Introducing 
a notation of a random variable x(t) standing for both of x p (t) and x u (t), we require the 
same form of stochastic differential equations for this x(t) as ([!]) and (fj) all the time, 

dx(t) = b(x(t),t)dt + dw(t) forward in time , (62) 
dx(t) = b*(x(t), t)dt + dw*(t) backward in time, (63) 

with the same properties for dw(t) as in([5]) and so on. Each sample path is described by 
x(t) as a whole, but has to be classified into either x p (t) or x u (t) at each t. Typically a 
sample path changes as, for example, x p {ti) — > x u {t2) — > • ■ ■ — > x p (t n ) as a result of repeated 
jumping processes. A sample path is said to be physically relevant at t if the sample is 
represented by x p (t), while it is not so if it is represented by x u (t). In other words, the 
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physical average at t is given by the average over ensemble of not all sample paths but only 
physically relevant sample paths at t. The notation <C f(x(t)) ^{x p (t)} is introduced to 
represent this conditional average for f(x(t)). Similarly the notations of other conditional 
averages such as < f(x(t)) »{ x „(t)} and < f(x(t)) »{x p (t)}u{x„(t)} are clear, in particular 

< f(x(t)) ^ {Xp(t)}u{Xu(t)} =< f( x (t)) > . (64) 

Again conditional averages related to many times can be introduced, e.g., <C 
f(x(t)) ^> { Xp (t+dt)}n{x u (t)} is supposed to represent the average over all the sample paths 
which are described by x u at t and x p at t + dt. 

Let us summarize the "dynamical" rule for stochastic jumping processes between p and 
u. The rules are given as follows (Fig. 2): (i) For forward time direction, in case of U < 0, 



x p (t) 



x u (t + dt) with the probability of —2U(x p (t),t)/hdt, 

x p (t + dt) with the probability of 1 + 2U(x p (t),t)/hdt 
Xu{t) — > x u(t + dt) with the probability of 1, (65) 



and in case of U > 0, 
x u (t) — ► < 



x p (t + dt) with the probability of 2U(x u (t),t)/hdt, 

x u (t + dt) with the probability of 1 — 2U(x u (t),t)/Kdt 
Xp{t) — > x p (t + dt) with the probability of 1. (66) 



(ii) For backward time direction, in case of U < 




x p (t) 
x u (t) 



Xp (t 


- dt) 


with the probability of —2U(x u (t),t)/hdt, 




x u (t 


- dt) 


with the probability of 1 + 2U(x u (t),t)/hdt 




Xp (t 


-dt) 


with the probability of 1, 


(67) 


Xu (t 


- dt) 


with the probability of 2U(x p (t),t)/hdt, 




Xp (t 


-dt) 


with the probability of 1 — 2U(x p (t),t)/hdt 




x u (t 


-dt) 


with the probability of 1. 


(68) 
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Note that for forward time direction the jumping process from x p to x u is allowed and the 
reverse process is forbidden where U < 0, and vice versa where U > 0, and that when U is 
non-positive everywhere, the number of sample paths described x p (t) decreases and that in 
x u (t) increases as t goes, the total number being conserved. Regardless of the indices of p 
and u, each sample path is a stochastic process described by fltx2| ) (or (]S3D) 

To prove the equivalence between the Fokker-Planck equation (|58|) and the stochastic 
differential equation (|62|) with the jumping rules (|65|) and (|66|), we calculate, for example, 
d < /(&(<)) 



dt 



1 r 

dt 
1 

dt 



< /(x(t + dt)) -> {Xp[t+ dt)} ~ < /(*(*)) >{*„(*)} 

< /(x(t + dt)) - f(x{t)) -> {Xpit +dt)}n{x p (t)} 
+ < /(z(t + dt)) ^>{ Xp (t+dt)}n{x u (t)} 

- < f(x(t)) >{xu(t+di)}n{x p (t)} 



(69) 



with 



< /(z(t + dt)) - /(x(t)) > {xp (i+d i )}n{ a;p (t)} 
= cit /ds /(x) 



9 h d 2 \ 

—b(x,t) + — — P(x,t)+o(dt 2 ), 
ox 2m ox z / 



(70) 



2U(x t) 

< f(x(t + dt)) > {Xp{t+d t)}n{x u (t)}= dt jdxfix) y ' P(x,t)e(U(x,t)) + o(dt 2 ), (71) 



and 



< /(x(t)) > { , u (^) } n { .. pW} = / dx / (g) 2[7(X ' l) P(x, t)0(-U(x, t)) + o(dt 2 ) . (72) 



These equations fl69D^fl72|) follows (0). The equivalence between the Fokker-Planck ap- 
proach and the approach of the stochastic differential equation (|62| ) with the stochastic 
jumping process (^) and ( |66| ) has been shown for forward direction. Similarly the equiva- 
lence between the two approaches can be proven for backward time direction. 

As for the dynamical condition, we do not modify the original Nelson's formulation. 
When the mean time derivatives Df(t) and _D*/(t) are concerned, there may be some am- 
biguity with respect to taking expectation. Here we will follow the argument given above 
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((49p. We define the "mean balanced acceleration" through the "mean time derivatives" as 
usual but for the stochastic process without any jumping process. We simply consider a 



stochastic process governed by (|62[) and (|63| ) all time. This leads to the "Newton-Nelson 
equation" in ( |T6"D in the present case. 

The combination of the equations (^T|)+i(|I^) leads to 

ox mip' ot 2 m ip' ox 1 m 

where the relation fll8|) is used. Again the relation between ip' and the solution of (|57|) ip is 
given as 

^ = ^'exp(-— J r)(s)ds), (74) 



and 



h d 

b(x,t) = — (Im + Re)— lnip(x,t), (75) 
m ox 

h d 

bJx,t) = — (Im-Re)— lnip(x,t), (76) 
m ox 

P(x,t) = mx } t)\ 2 . (77) 



are satisfied. 



V. NUMERICAL ANALYSIS 



Now we can perform the numerical analysis of the effects of the optical potential and 
channel coupling on the tunneling time, using above generalized Nelson's approach. 
First, we discuss one-dimensional system with a static square well optical potential, 



in I (x< 0), 

V{x) = \v -iU in II (0 < x < d), 
in III (d < x) 

(Fig. 3). We set the solution of the Shroodinger equation 



(78) 
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ih—ip{x,t) 



2m dx 2 



+ V(x)]ip(x,t) 



(79) 



as 



i>(x,t)= A(k)ip k (x)e' l ^dk (80) 

J— oo 

with a coefficient function A(k) and E = It is well-known that (pk(x) is written as 

in I, 



where 



C k e KX + D k e~ KX 



ikx 



in II, 



in III, 



/2m(V -iU -E) 

«=- Z =K R -iK! (K/>0), 

ft 



(81) 



(82) 



/2m(V - itTb - #o) 

( /.',, = = K m - IKiq [K I0 > 0), ) 



and R k i T k , C k and D k are given as 



Rk 




—i(n 2 + k 2 )smh.nd 


n 




2kKe~ ikd 




= B 




Ck 




k(n + ik)e~ Kd 


D k 




k(n - ik)e Kd 



(83) 



(84) 



B 



1 



2kKcosh.Kd + i{k 2 — k 2 )sinh.Kd 
We take a Gaussian form with its center at k = ko, or 

f (k -k) 2 ) 



A(k) = A ko (k) = Cexp 



4a 2 



with a normalization constant C. We put here a = and Vb = 5i?o 
solution, we calculate numerically (|62|), ( |63"D and (^) ~ (^) . 



(85) 



(86) 



;|^) 2 . Using this 
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Figure 4 shows the three typical sample paths calculated by ( p2|) and ( p5|) and (p6|). 
There is a sample path x(t) which change it's property from "physical" to "unphysical" in 
the tunnel region. 

Figures 5 and 6 shows the parameter ^ v.s. the average of passing time t p , calculated 
by (0), fl6?D and (|68|) . See the details of this "backward time evolution method" in our 
previous work p0|| . Generally, r p decrease as the ^ become larger. Let us estimate r p 
analytically on the W.K.B. like approximation. If we can write the wave function in the 
tunnel region II as 

ip{x,t) ~ tp(x) ~ C'exp(-K x) = exp{-(K jR0 - iK I0 )x}, (87) 
the drift of (|63|) becomes 



0* = — (K/o + Kflo) (H hO )), K =- 7 • (88) 

777 777 K R0 K R0 ft 

from (|76|) . In these cases, the "backward" time evolution of the distribution function Pt{x, t), 
which has an "initial" distribution 5(x — d), is written as 

dP T (x,t) r % . k io . d h d 2 2 tt1t ^ . . 

and we can get the solution of (|89|) easily 

^0(1+-^-) 



M*,t) = \I—*M- r 1 ^ ( t<0 )- (9°) 



2m7r Ax-d + ^t) 2 2U 

d 2 



-ht " L ft 

2m 



There are two characteristic time intervals in this solution. One is the diffusion time td ~ ^\ 
for which the distribution sizes up to the potential width d. The other is the current time 

md md k io 

^fii^TTsTfii^ 1 -^' (91) 

for which the peak of the distribution moves from x — d to x — 0. Of course, the approxi- 
mation of (^) is justified when Kod is much larger than 1, and this leads us to the relation 
of 

U > t c , (92) 
19 



and the time interval t c becomes the passing time in this extreme case. Note that this t c 
has tendency of decreasing as the ^ become larger. 

Second, we discuss a one-dimensional system with a static square well potential and 
2-cannel coupling, or the case of the Shrodinger equation for this problem written down as 



d 


in 






4>2 





2m dx 2 

u 



u 



2m dx 2 



+ v 



V>2 



(93) 



V and U are supposed to 



and 










in I 


< o), 




V(x) 


= < 


V 


in II 


(0 < x < d), 


(94) 









in III 


(d < x), 








f 




in I 


(x < 0), 




U(x) 




U 


in II 


(0 < x < d), 


(95) 









in III 


(d < x). 





Figure 7 shows the schematical illustration of our simulation. We can diagonalize 





where 



d 


1>+ 




l dt 










i 






71 







as 



(96) 



(ipi + fa) and ip. 



V2 



(97) 



and we can write down the time-dependent solution of ip + and ip- easily as same as (|30"D or 



ip±(x,t) 



A(k)^ ±k {x)e- l ^dk 



(98) 



with a Gaussian coefficient function A(k) and E = hr~. <f±k{ x ) is (HD substituted k for 



K± 



2m{V ±U -E) 
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(99) 



Figure 8 shows the some typical sample paths calculated by (p4|) , (|37|) and (^) . There is 
a path which changes its index from 1 to 2 in the passage through the tunneling region. t i=12 
Figs. 9 and 10 are the averages of the passing times over the sample paths which belong to 
{xi(t)} at t — > oo. We can see in Fig. 10 that there is a critical parameter value of Vo ^° = 1 
in the behavior of t\ and t^. This is understood as following; In the case of Vo ^ > 1, the 
"— " channel, which is dominant in the tunnel region II in comparison with"+" one, is not 
tunneling channel and it describes a particle which goes over the potential. Regardless of 
xi(t) and X2(t), the time spent in the "potential region" is expected to be agree with the 
one which is expected from classical mechanics or 

md 



where k_ 



2m(E -V + U ) 



(100) 



hk_ h 

This is seen in the Fig. 9 too. On the other hand, in the case of v °^° < 1, we can 
approximate the wave functions ip + and ip- in the thick tunnel region to 

(101) 














C- exp(— K_ x) 



and ipi and to 



where 





1 




exp(— k_ x) 


^2 




-c. 


_ exp(— k_ x) 


K 




-Uq- E ) 


-o — 




h 



(102) 



(103) 



So we can estimate "passing time" of both channels (1 and 2) at 

md 

TlK_n 



(104) 



likewise ([FT]). 
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VI. SUMMARY AND COMMENTS 



In this paper, we have analyzed the effects of inelastic scattering on the tunneling time 
theoretically, using generalized Nelson's quantum mechanics. This generalization enable us 
to describe quantum system with optical potential and channel couplings in a real time 
stochastic approach. In this formalism, the space-time development of dynamical variable, 
e.g. coordinate of particle, is described by a definite path determined stochastically. Each 
sample path has a definite form of trajectory in space-time diagram, while a physical quantity 
averaged over the ensemble of these sample paths recovers the effect of quantum coherence. 
This is true even in the Young's double slits interference experiment. Nelson's quantum 
mechanics gives each definite trajectory and the ensemble of it, but it does not predict 
which path is selected when one wants to measure the position of a particle. In this sense, 
this "real time stochastic process approach" seems to give us a new insight into quantum 
mechanics beyond Copenhagen interpretation. On the other hand, the effects of more general 
cases (many body system, environment, temperature and so on) are subjects in a future and 
this work would be the first step to that study. 

Recent experimental data of tunneling time using the neutron spin echo shift through 
the magnetic films ]T5[ seem to agree with the simulation based on our approach |21| and 



this study will be reported in near future. 
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Figure caption 



FIG. 1. Schematical illustration of the "dynamical" rule for stochastic jamping process 
between two cannels. 

FIG. 2. Schematical illustration of "dynamical" rule for stochastic jamping process between 
physical and unphysical sector. 

FIG. 3. Schematical illustration of one-dimensional optical barrier tunneling. 

FIG. 4. The three typical sample paths in the optical potential case. 

FIG. 5. The mean value of r p versus j& (thin potential cases). 

FIG. 6. The mean value of r v versus ^ (thick potential cases). 

FIG. 7. Schematical illustration of one-dimensional scattering with channel coupling. 

FIG. 8. The three typical sample paths with the the channel coupling. 

FIG. 9. The mean values of t\ and t 2 versus ^ (thin potential cases). 

FIG. 10. The mean values of t\ and t 2 versus ^ (thick potential cases). 
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